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Abstract In an analysis of XMM-Newton archival observations of the bright Local Group spiral galaxy M 31 we study the 
population of X-ray sources (X-ray binaries, supernova remnants) down to a 0.2^1.5 keV luminosity of 4.4x 10 34 erg s" 1 . EPIC 
hardness ratios and optical and radio information are used to distinguish between different source classes. The survey detects 
in an area of 1.24 square degree 856 sources. We correlate our sources with earlier M 31 X-ray catalogues and use information 
from optical, infra-red and radio wavelengths. As M 3 1 sources we detect 2 1 supernova remnants (SNR) and 23 SNR candidates, 
18 supersoft source (SSS) candidates, 7 X-ray binaries (XRBs) and 9 XRB candidates, as well as 27 globular cluster sources 
(G1C) and 10 G1C candidates, which most likely are low mass XRBs within the G1C. Comparison to earlier X-ray surveys 
reveal transients not detected with XMM-Newton, which add to the number of M 31 XRBs. There are 567 sources classified 
as hard, which may either be XRBs or Crab-like SNRs in M 31 or background AGN. The number of 44 SNRs and candidates 
more than doubles the X-ray detected SNRs. 22 sources are new SNR candidates in M 31 based on X-ray selection criteria. 
Another SNR candidate may be the first plerion detected outside the Galaxy and the Magellanic Clouds. On the other hand, 
six sources are foreground stars and 90 foreground star candidates, one is a BL Lac type active galactic nucleus (AGN) and 36 
are AGN candidates, one source coincides with the Local Group galaxy M 32, one with a background galaxy cluster (GC1) and 
another is a GC1 candidate, all sources not connected to M 3 1 . 
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1. Introduction 

In the XMM-Newton survey of the Local Group Sc galaxy 
M 33, Pietsch et al. (2004, hereafter Paper I) detected 408 
sources in a 0.8 square degree field combining the counts of all 
EPIC instruments. The use of X-ray colours and optical and ra- 
dio information allowed them to identify and classify the X-ray 
sources and proved to be efficient in separating super-soft X- 
ray sources (SSSs) and thermal supernova remnants (SNRs) in 
M 33 from Galactic stars in the foreground and "hard" sources. 
These hard sources may be either X-ray binaries (XRBs) or 
Crab-like SNRs in M 33 or active galactic nuclei (AGN) in the 
background of the galaxy. 

The Andromeda galaxy M 31 is located at a similar dis- 
tance as M 33 (780 kpc, Holland 1998; Stanek & Garnavich 
1998, i.e. 1" corresponds to 3.8 pc and the flux to luminos- 
ity conversion factor is 7.3xl0 49 cm 2 ) and - compared to the 
near face-on view of M 33 - is seen under a higher inclina- 
tion (78°). The optical extent of the massive SA(s)b galaxy can 
be approximated by an inclination-corrected D25 ellipse with 



* XMM-Newton is an ESA Science Mission with instruments and 
contributions directly funded by ESA Member States and the USA 
(NASA). 

** Table 2 is available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u- 
strasbg.fr/cgi-bin/qcat ?J/A+ A/ 



large diameter of 153'3 and axes ratio of 3.09 (de Vaucouleurs 
et al. 1991; Tully 1988). With its moderate Galactic foreground 
absorption (N H = 7xl0 20 cm" 2 , Stark et al. 1992) M 31 is well 
suited to study the X-ray source population and diffuse emis- 
sion in a nearby spiral similar to the Milky Way. As M 31 is 
seen through a comparable Galactic absorbing column as M 33 
one can use the same methods and similar source selection cri- 
teria that proved to be successful in the M 33 analysis. On the 
other hand, there is a big number of low mass XRBs (LMXBs) 
identified as bright X-ray sources in M 3 1 from earlier obser- 
vations (see below). The properties of these sources may help 
to better classify some of the "hard" sources in M 33. 

M 3 1 was a target for many previous imaging X-ray mis- 
sions. The Einstein X-ray observatory detected 108 individual 
X-ray sources brighter than ~5xl0 36 erg s" 1 , 16 of which were 
found to vary between Einstein observations (van Speybroeck 
et al. 1979; Collura et al. 1990; Trinchieri & Fabbiano 1991). 
The sources were identified with young stellar associations, 
globular clusters (i.e. LMXBs) and SNRs (see e.g. Blair et al. 
1981; Crampton et al. 1984). With the ROSAT HRJ, Primini 
et al. (1993) reported 86 sources brighter than ~10 36 erg s _1 in 
the central area of M 31, nearly half of which were found to 
vary when compared to previous Einstein observations. With 
a separation of about one year between the M 3 1 surveys, the 
ROSAT PSPC covered the entire galaxy twice and detected al- 
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together 560 X-ray sources down to a limit of ~5xl0 erg s 
(Supper et al. 1997, 2001, hereafter SHP97, SHL2001). The 
intensity of 34 of the sources varied significantly between the 
observations, and SSSs were established as a new class of M 31 
X-ray sources (see also Kahabka 1999). 

With the new generation of X-ray observatories, Chandra 
and XMM-Newton, up to now only parts of M 31 were sur- 
veyed. Deep Chandra ACIS-I and HRC observations of the 
central region (covered areas of 0.08 and 0.27 square degree) 
resolved 204 and 142 X-ray sources, respectively (Garcia et al. 
2000; Kong et al. 2002b; Kaaret 2002). Three M 31 disk fields, 
spanning a range of stellar populations, were covered with short 
Chandra ACIS-S observations to compare their point source 
luminosity functions to that of the galaxy's bulge (Kong et al. 
2003a). A synoptic study of M 31 with the Chandra HRC cov- 
ered in 17 epochs "most" of the disk (0.9 square degree) in 
short observations and resulted in mean fluxes and long-term 
light curves for the 166 objects detected (Williams et al. 2004). 
In these observations, several M 31 SNRs were spatially re- 
solved (Kong et al. 2002a, 2003b) and bright XRBs in globular 
clusters and SSSs and quasisoft sources (QSSs) could be char- 
acterized (Di Stefano et al. 2002, 2004; Greiner et al. 2004). 
During the XMM-Newton guaranteed time program there were 
four observations of the central area of M 3 1 and three aimed 
at the northern and two at the southern disk. These observa- 
tions were used to investigate the bright and variable sources 
and diffuse emission (Shirey et al. 2001; Osborne et al. 2001; 
Trudolyubov et al. 2001, 2004), and to derive source luminosity 
distributions (Trudolyubov et al. 2002a). In addition, the time 
variability and spectra of several individual XRBs have been 
studied in detail (e.g. Trudolyubov et al. 2002b; Barnard et al. 
2003a,b; Mangano et al. 2004). 

Here we present X-ray images and a source catalogue for 
the archival observations of M 3 1 combining the three XMM- 
Newton EPIC instruments and using only times of low back- 
ground. For the source catalogue and source population study 
of these observations we analyzed the individual pointings and 
the merged data of the central area simultaneously in five en- 
ergy bands in a similar way as described in our M 33 analysis 
in Paper I. The covered area of 1.24 square degree and lim- 
iting sensitivity is a significant improvement compared to the 
Chandra surveys, however, up to now only covers about 2/3 of 
the optical M 3 1 extent (D25 ellipse) with a rather inhomoge- 
neous exposure. 

2. Observations and data analysis 

Table 1 summarizes the archival XMM-Newton (Jansen et al. 
2001) EPIC (Struder et al. 2001; Turner et al. 2001) observa- 
tions which at least partly overlap with the inclination corrected 
optical M 3 1 D25 ellipse. For each observation we give the field 
name and the abbreviation used in the text for the M 31 ob- 
servation (Col. 1), the observation identification (2), date (3), 
pointing direction (4,5) and systematic offset (6), as well as fil- 
ter and exposure time after screening for high background for 
EPIC PN (7,8), MOS1 (9,10), and MOS2 (11,12). For creat- 
ing the full field colour image we had to apply an even more 



stringent background screening which resulted in the shorter 
exposure times given in brackets. 

In the XMM-Newton observations the EPIC PN and MOS 
instruments were operated in the full frame mode providing a 
time resolution of 73.4 ms and 2.6 s, respectively. The medium 
filter was in front of the EPIC cameras in all but the obser- 
vations c4, si and s2 which were performed with the thin fil- 
ter. For creating the merged images and source detection we 
used medium and thin filter observations together. This proce- 
dure can be justified as the difference in absorption between the 
medium and soft filters corresponds to an 7Y H = l.lxlO 20 crrT 2 , 
significantly lower than the Galactic foreground absorption. 

We carefully screened the event files for bad CCD pixels 
remaining after the standard processing. To create a homoge- 
neous combined colour image with similar background level 
for all fields we had to carefully screen the data for times 
of high background using the high energy (7-15 keV) back- 
ground light curves provided by the SAS tasks epchain and 
emchain. For the images that were used for the source de- 
tection procedures in the individual fields we adopted a less 
stringent background screening which due to the then longer 
exposures allowed us to detect fainter sources. The good time 
intervals (GTI) were determined from the higher statistic PN 
light curves and also used for the MOS cameras. Outside the 
PN time coverage, GTIs were determined from the combined 
MOS light curves. The corresponding low background times 
for the individual observations are listed for the PN and MOS 
cameras in Table 1 . 

The archive contains four observations of the centre area 
of M 3 1 separated by half a year adding up to low background 
exposures for source detection of 124.6, 133.1 and 133.2 ks for 
EPIC PN, MOS1 and MOS2 cameras, respectively. In addition 
there are two pointings in the southern and three in the north- 
ern disk of M 31 that aimed for 60 ks exposure. Due to high 
background, for two of the observations the time usable for de- 
tection is significantly shorter (s2, n2). In addition the archive 
contains four shorter pointings into the halo of M 31 (hl-h4, 
PI Di Stefano). These observations mainly contain foreground 
and background objects as they were pointed far of the disk of 
the galaxy. Only observation h4 (1 1 ks low background) partly 
covers the disk and therefore is included in our analysis. In to- 
tal, the observations in our analysis cover an area of 1 .24 square 
degree. 

The data are treated in a similar way as the M 33 data 
described in Paper I. We used five energy bands: (0.2-0.5) 
keV, (0.5-1.0) keV, (1.0-2.0) keV, (2.0-4.5) keV, and (4.5-12) 
keV as band 1 to 5. We intentionally split the (0.5-2.0) keV 
band used in the 1XMM XMM-Newton Serendipitous Source 
Catalogue 1 to get on average a more homogeneous distribution 
of the source counts to the energy bands which leads to a better 
spread of the hardness ratio values and allows effective source 
classification (see Paper I and Sect. 5). 

For PN we selected only "singles" (PATTERN = 0) in band 
1, for the other bands "singles and doubles" (PATTERN < 4). 
For MOS we used "singles" to "quadruples" (PATTERN < 12). 

1 Prepared by the XMM-Newton Survey Science Centre 
Consortium (http : //xmmssc-www. star . le . ac . uk/) 
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Table 1. XMM-Newton log of archival M 31 observation overlapping with the optical D 2 5 ellipse (proposal numbers 010927, 011257 and 015158). 
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M31 field 


Obs. id. 


Obs. dates 


Pointing 


direction 


Offset * 


EPIC PN 


EPICMOS1 


EPIC MOS2 








RA/DEC (J2000) 




Filter* 




Filter* 


Tip 


Filter* 


T f exp 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Centre 1 (cl) 


0112570401 


2000-06-25 


0:42:36.2 


41:16:58 


-1.9, +0.1 


medium 


26.40(23.60) 


medium 


29.92(29.64) 


medium 


29.91(29.64) 


Centre 2 (c2) 


0112570601 


2000-12-28 


0:42:49.8 


41:14:37 


-2.1, +0.2 


medium 


9.81( 5.85) 


medium 


12.24( 6.42) 


medium 


12.24( 6.42) 


Centre 3 (c3) 


0112570701 


2001-06-29 


0:42:36.3 


41:16:54 


-3.2,-1.7 


medium 


27.65(21.83) 


medium 


27.65(23.85) 


medium 


27.65(23.86) 


North 1 (nl) 


0109270701 


2002-01-05 


0:44:08.2 


41:34:56 


-0.3, +0.7 


medium 


54.78(48.57) 


medium 


57.31(55.68) 


medium 


57.30(55.67) 


Centre 4 (c4) 


0112570101 


2002-01-06/07 


0:42:50.4 


41:14:46 


-1.0, -0.8 


thin 


60.79(48.11) 


thin 


63.31(52.87) 


thin 


63.32(52.86) 


South 1 (si) 


0112570201 


2002-01-12/13 


0:41:32.7 


40:54:38 


-2.1,-1.7 


thin 


53.45(47.00) 


thin 


53.76(51.83) 


thin 


53.73(51.84) 


South 2 (s2) 


0112570301 


2002-01-24/25 


0:40:05.0 


40:34:38 


-1.1,-0.3 


thin 


29.31(22.23) 


thin 


29.37(24.23) 


thin 


29.38(24.24) 


North 2 (n2) 


0109270301 


2002-01-26/27 


0:45:28.6 


41:55:26 


-0.3,-1.5 


medium 


26.27(22.86) 


medium 


26.54(25.22) 


medium 


26.55(25.28) 


North 3 (n3) 


0109270401 


2002-06-29/30 


0:46:31.1 


42:17:32 


-2.3,-1.7 


medium 


51.39(41.28) 


medium 


59.69(45.39) 


medium 


59.69(45.53) 


Halo 4 (h4) 


0151580401 


2002-02-06 


0:46:15.9 


41:20:25 


+0.3, +0.0 


medium 


11.29(11.29) 


medium 


12.91(12.91) 


medium 


12.92(12.92) 



Notes: 

* : Systematic offset in Ra and DEC in arcsec determined from correlations with 2MASS, USNO-B1 and Chandra catalogues 
+ : all observations in full frame imaging mode 

f : Exposure time in units of ks after screening for high background used for detection, for colour image in brackets (see text) 



To avoid background variability over the PN images we omit- 
ted the energy range (7.2-9.2) keV in band 5 where strong flu- 
orescence lines cause higher background in the outer detector 
area (Freyberg et al. 2004). To convert source count rates in 
the individual bands to fluxes we used count rate to energy 
conversion factors (ECF) for PN and MOS observations calcu- 
lated for the epoch of observation and filters used, assuming the 
same spectrum as for the first XMM-Newton source catalogue, 
i.e. a power law spectrum with photon index 1.7 absorbed by 
the Galactic foreground column of 7xl0 20 cirT 2 . As has been 
demonstrated in Table 2 of Paper I, ECF values only vary by 
about 20% for PN in energy bands 1 to 4 for different spectra 
like an absorbed 1 keV thin thermal model typical for SNRs or 
a 30 eV black body model typical for SSSs. The same is true for 
EPIC MOS, with the exception of the black body ECF in the 
0.2-0.5 keV band where the MOS sensitivity is much lower. 
For the detection in the merged centre images (combining ob- 
servations with thin and medium filter) we used averaged ECFs 
which may lead to wrong flux estimates of at maximum 5%. 

For most sources band 5 just adds noise to the total count 
rate. If converted to fluxes due to the high ECF this noise often 
dominates the total flux. To avoid this problem we calculated 
count rates and fluxes for detected sources in the "XID" (0.2- 
4.5) keV band (bands 1 to 4 combined). While for most sources 
this is a good solution for extremely hard or soft sources there 
may still be bands just adding noise. This then leads to rate 
and flux errors that seem to wrongly indicate a lower source 
significance. A similar effect occurs for the all instrument rates 
and fluxes if a source is mainly detected in one instrument (e.g. 
soft sources in PN). 

To classify the source spectra we computed four hardness 
ratios from the source count rates. These hardness ratios and 
errors are defined as 

HRi = M and EHRi = 2 ^pgp^ , 

for ; = 1 to 4, where B, and EBi denotes count rates and 
corresponding errors in band i as defined above. In the standard 
source detection products hardness ratios and count rates in in- 



dividual energy bands are not combined for all instruments. To 
improve statistics we calculated for all energy bands "all EPIC" 
count rates, fluxes and hardness ratios. We are aware that these 
products have some jitter in their meaning, as the relative in- 
tegration times for individual sources in the EPIC instruments 
differ. Due to the better signal in the combined products they 
are still very valuable for source classification. For investiga- 
tions needing accurate calibration the products for the individ- 
ual instruments should be used that are also given. 

We created for PN, MOS1 and MOS2 in each of the 5 en- 
ergy bands mentioned above images, background images, ex- 
posure maps (without and with vignetting correction), masked 
them for acceptable detector area. For PN the background 
maps contain the contribution from the "out of time (OOT)" 
events (parameter withootset=true in task esplinemap). 
In contrast to our M 33 raster, the M 31 images (with the 
exception of the four observations towards the M 31 cen- 
tre, cl, c2, c3, c4) only overlap at the edge of the field of 
view (FOV). We therefore searched for sources in each of 
these fields individually and only merged the centre fields 
for source detection. For visualization purposes we created 
images of all observations individually, for the centre and 
for the full area merging all EPIC instruments (see Sect. 4). 
To allow an easy merging of the galaxy centre images, we 
re-calculated for the events of all contributing observations 
projected sky coordinates with respect to reference position 
RA=00 h 42 m 45 s .0, DEC=+41°15'00"(J2000), for the full field 
images with respect to reference position RA=00 h 43 m 24:0, 
DEC=+41°24'00"(J2000), respectively. 

To create event lists and images we used the lat- 
est calibration products for the linearisation of the EPIC 
MOS CCDs (MOS*_LINCOORD_0017.CCF) and EPIC bore- 
sight (XMM.BOPvESIGHT.0018.CCF). We verified with the 
source-rich M 31 centre pointings that PN and MOS source 
positions coincide to better than 0.5". On the other hand there 
were still significant offsets between the observations that had 
to be corrected for before merging. For the centre observa- 
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tions these offsets were determined from source lists of the in- 
dividual observations. We used the USNO-B1, 2MASS, and 
Chandra catalogues to define an absolute reference frame. The 
offets applied are listed in column 6 of Table 1 . This finally re- 
sulted in a residual systematic position error of less than 0.5". 

The images for PN, MOS 1 and MOS2 do not fully overlap 
at the borders due to the different FOV of the EPIC instruments. 
Nevertheless, we used the full area for source detection. 

The data analysis was performed using tools in the SAS 
v6.0 and some later versions from the development area as spe- 
cially mentioned, EXSAS/MIDAS 03OCT_EXP, and FTOOLS 
v5.2 software packages, the imaging application DS9 v3.0b6 
and spectral analysis software XSPEC vl 1.2. 

3. Source catalogue 

We searched for sources using simultaneously 5x3 images 
(5 energy bands and PN, MOS1 and MOS2 camera). For the 
pointings into the disk and halo of M 31 this procedure was 
applied to the individual observations. The centre pointings 
strongly overlap and we therefore merged the images to reach 
higher detection sensitivity. A preliminary source list created 
with the task boxdetect with a low likelihood threshold was 
used as starting point for the task eraldetect v. 4.40. To re- 
solve sources that overlap due to the PSF, we used parameter 
multisourcefit=4. We selected sources equivalent to a sin- 
gle camera likelihood of 10 by accepting sources in the com- 
bined fit which have a likelihood above 7. This procedure pro- 
duced acceptable results when checked against the smoothed 
images (see Sect. 4). Two emission peaks were resolved into 
several sources with no sign of multiplicity from the smoothed 
images. In a specific detection run we found these sources (747 
and 832) with a high likelihood for extent and that they can 
be characterized by a King beta model with a core radius of 
3071 ±073 and 2873+172, respectively. We rejected spurious 
detections in extended emission around the centre. Finally we 
merged the source lists of the centre and individual pointings 
to the disk and ordered the sources according to increasing 
right ascension. For sources detected in more than one observa- 
tion we included the parameters of that source in the catalogue 
that was found with the higher detection likelihood. A recur- 
rent transient source was reported from Chandra and XMM- 
Newton observations (Garcia et al. 2000; Osborne et al. 2001; 
Williams et al. 2004). It is only visible during the short obser- 
vation c2 and was not significantly detected above the central 
diffuse emission by the detection procedure in the merged cen- 
tre observation. We therefore derived the source parameters of 
this source from just observation c2, and added them at the end 
of the catalogue. In total we detect 856 sources in the field. 

The source parameters are summarized in Table 2 (EPIC 
combined products and products for EPIC PN, MOS1 and 
MOS2, separately) which is available in electronic form at the 
CDS. 

Table ?? gives the source number (Col. 1), detection field, 
from which the source was entered into the catalogue (2), 
source position (3 to 9) with lcr uncertainty radius (10), likeli- 
hood of existence (11), integrated PN, MOS1 and MOS2 count 
rate and error (12,13) and flux and error (14,15) in the (0.2- 



4.5) keV XID band, and hardness ratios and errors (16-23). 
Hardness ratios are calculated only for sources for which at 
least one of the two band count rates has a significance greater 
than 2cr. Errors are the properly combined statistical errors in 
each band and can extend beyond the range of allowed values 
of hardness ratios as defined previously (-1 .0 to 1 .0). The EPIC 
instruments contributing to the source detection, are indicated 
in the "Val" parameter (Col. 24, first character for PN, second 
MOS1, third MOS2) as "T", if inside the field of view (FOV), 
or "F", if outside of FOV. There are 52 sources at the periphery 
of the FOV where only part of the EPIC instruments contribute. 
The positional error (10) does not include intrinsic systematic 
errors which amount to 075 (see above) and should be quadrat - 
ically added to the statistical errors. 

Table 2 then gives for EPIC PN exposure (25), source ex- 
istence likelihood (26), count rate and error (27,28) and flux 
and error (29,30) in the (0.2-4.5) keV XID band, and hard- 
ness ratios and error (31-38). Columns 39 to 52 and 53 to 66 
give the same information corresponding to Cols. 25 to 38, but 
now for the EPIC MOS1 and MOS2 instruments. Hardness ra- 
tios for the individual instruments were again screened as de- 
scribed above. From the comparison of the hardness ratios de- 
rived from integrated PN, MOS 1 and MOS2 count rates (Cols. 
16-23) and the hardness ratios of the individual instruments 
(Cols. 31-38, 45-52 and 59-66) it is clear that combining the 
instrument count rate information yielded significantly more 
hardness ratios above the chosen significance threshold. 

Column 67 shows cross correlations with M 3 1 X-ray cata- 
logues in the literature (see detailed discussion in Sect. 5). Only 
364 sources in our catalogue correlate with previously reported 
M 3 1 X-ray sources, i.e. our catalogue doubles the number of 
known X-ray sources in M 3 1 . 

In the remaining columns of Table 2, we give cross cor- 
relation information with sources in other wavelength ranges, 
which is further described in Sect. 6. We only want to men- 
tion here, that we used the foreground stars, globular cluster 
sources and candidates, to verify the assumed source position 
errors. All but six of the 33 identifications and 100 candidates 
are located within the 3cr statistical plus systematic positional 
error given above. 

The faintest sources detected have an XID band flux of 
6.0xl0~ 16 erg cm" 2 s" 1 . The brightest source (50, XID band 
flux of 4.5xl0~ 12 erg cirT 2 s _1 ) is identified as AGN in the 
M 31 field. The brightest M 31 source (291, an XRB) has an 
XID band flux of 3.9xl0~ 12 erg cirT 2 s" 1 . The XMM-Newton 
detected sources in M 3 1 therefore cover an absorbed luminos- 
ity range of 4.4xl0 34 erg s -1 to 2.8xl0 38 erg s _1 in the XID 
band. 

4. Images 

For Figs. 1, 3, 4, we smoothed the images and exposure maps 
used for detection with a Gaussian (FWHM of 5", 10" and 
20", respectively). For PN we subtracted OOT images, that 
were masked and smoothed in the same way as the images. 
For the EPIC combined images we added the images of the 
individual cameras scaled according to the background in the 
individual energy bands. For the colour image Fig. 4, we added 
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Figure 1. Logarithmically-scaled XMM-Newton EPIC low background image with a pixel size of 1 arcsec 2 of the M 31 centre observations 
combining PN and MOS1 and MOS2 cameras in the (0.2-4.5) keV XID band. The data are smoothed with a Gaussian of FWHM 5" which 
corresponds to the point spread function in the centre area of the combined observations. The image is corrected for un-vignetted exposure and 
masked for exposures above 5 ks for the individual cameras. Contours are at (1,2,4,8, 16,32) x 1CT 6 ct s -1 pix~' including a factor of two 
smoothing. Sources from the catalogue are marked in the outer area. The inner area is shown in detail in Fig. 2. 



images for the individual bands as needed, masked them to a 
total un-vignetted exposure of 5 ks and indicated the optical 
extent of M 3 1 by the D25 ellipse. 

The (0.2-^.5) keV XID band images (Figs. 1, 3) give an 
overview on the sources detected in the XMM-Newton analy- 
sis. To better visualize faint structures we added contours. By 
comparing images in the different energy bands it is clear that 
many sources only show up in some of the images, a fact that 



indicates spectral diversity and is further quantified in the dif- 
ferent hardness ratios of the sources. This fact can be visual- 
ized even more clearly in the combined EPIC colour image 
where we coded the (0.2-1.0) keV band in red, the (1.0-2.0) 
keV band in green and the (2.0-12) keV band in blue (see 
Fig. 4). The image is a demonstration of the colourful X-ray 
sky. SSSs, thermal SNRs and foreground stars appear red or 
yellow, XRBs, Crab-like SNRs and AGN green to blue. Bright 
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Figure 2. Inner area of Fig. 1. Contours are at (1, 2, 4, 8, 16, 32, 64, 128, 256)x 10 6 ct s 1 pix 1 . Left: Sources from the XMM-Newton catalogue 
are marked as 30"x30" squares. Right: Sources from the Chandra analysis of Kong et al. (2002b) and Williams et al. (2004) are marked. 



diffuse emission fills the bulge and fainter emission the area of 
the disk (yellow and red, see also Fig. 1). A detailed analysis of 
this emission is outside the scope of this paper. XMM-Newton 
analysis of the diffuse emission of the M 3 1 bulge region was 
reported by Shirey et al. (2001) and Takahashi et al. (2004) and 
of the northern disk by Trudolyubov et al. (2004). 

5. Cross-correlation with other M 31 X-ray 
catalogues 

In this section we discuss the cross-correlation of the XMM- 
Newton detected sources with sources reported in earlier X-ray 
catalogues. All correlations (together with other X-ray infor- 
mation like variability (v) or transient nature (t), reported in 
these catalogues as well as extent (ext) detected in this work) 
are indicated in the XID column of Table 2 (Col. 67). 

From the 108 Einstein sources reported in Trinchieri & 
Fabbiano (1991, hereafter TF91), 14 are outside the field 
covered by the XMM-Newton observations ([TF91] 6,9,11, 
13,15,38,81,86,95,97,98,102,106,108). From the remaining 94 
sources 22 are not detected in our analysis while they were 
clear Einstein detections ([TF91] 29,30,31,35,37,39,40,43, 
46,47,53,65,66,72,75,78,88,93,96,99,100,107) and have to be 
classified as transient. Three Einstein sources ([TF91] 10,70, 
103) are resolved by XMM-Newton into two sources. Several 
of the transient candidates ([TF91] 31,39,40,46,47,53,75,78, 
96) were already classified as transients or variable ([TF91] 37, 
88) after the ROS AT HRI and PSPC observations (Primini et al. 
1993, hereafter PFJ93, and SHL2001). On the other hand, the 
Einstein source [TF91] 84, classified as transient by SHL2001, 
was detected with XMM-Newton at an absorbed XID band lu- 
minosity of 1.5xl0 36 erg s" 1 , about a factor of 8 below the 
luminosity reported by TF91 and a factor of two below the 
SHL2001 upper limit. This source therefore either is a recur- 
rent transient or just highly variable. 



In the correlation with ROSAT detected sources we 
concentrate on the catalogue of the central HRI pointing 
(PFJ93) performed in July 1990, and the first (July 1991) 
and second (August 1992) PSPC survey of the full galaxy 
(SHP97,SHL2001). With the HRI 86 sources are detected 
within the central ~34' of M 31. The XMM-Newton observa- 
tions cover the full field and detect all but 7 of the sources. The 
missing sources ([PFJ93] 1,2,31,33,51,63,85) are not in con- 
fused regions and covered the luminosity range (0.2-2. 4)xl0 37 
ergs 1 . While [PFJ93] 1,2,85 were still detected in the first 
and/or second ROSAT PSPC survey, sources [PFJ93] 31,33,51, 
63 were no longer active. All ROSAT HRI sources should have 
clearly been detected as bright sources in the XMM-Newton 
survey and the missing sources therefore have to be classified 
as transients. 

The ROSAT PSPC surveys of M 31 contain 396 sources 
each, covering 6.3 and 10.7 deg 2 , an area much bigger than 
the 1.24 deg 2 of the XMM-Newton survey. In Table 3 we 
list sources of SHP97 and SHL2001 that are not detected by 
XMM-Newton. From ROSAT PSPC survey I and II, 156 and 
191 sources are outside the field covered by the XMM-Newton 
observations, 43 and 23 are not detected by XMM-Newton, 
respectively. For M 33 we found in Paper I that many of the 
sources with the lowest detection likelihoods LH were spurious 
detections. To search for similar effects in the M 31 ROSAT 
catalogues we arranged the ROSAT sources in the field, that 
are not detected by XMM-Newton according to LH. From the 
first ROSAT survey, we do not detect 12 out of 18 sources 
with LH below 12 and 9 out of 29 sources with LH below 15. 
These sources should have been detected by the deeper XMM- 
Newton observations if still at the ROSAT brightness. While in 
principle they all could have dimmed by such an amount that 
we do not detect them, this seems rather unlikely and proba- 
bly most of these sources were spurious detections. They may 



W. Pietsch et al.: XMM-Newton survey of M 31 







a 



10100 



□ 



E Do 

mn 



J3 



n m 

1 ™**, 



"7 



M 31 North 2 

XMVI-Newton 



| 
□ 



. -42100)00 »• 



□ S3 



13 



i 



us 

a 



a 



B « El* ... 

B 



3 sm im 



-4145.00 

0.2 -45 keV 
EPIC PM + M0S 



30 046100 30 15 O0 41 30 
-i—i i—i i—l i—i i—i i—l i—i i—i I I I i—i i—i l—i i—i i—i l—i i—i i— tJ 



TT3 

13 



I I 

M 31 Halo 4 

XMVI -Newton 











,9. n «3 



a* 

'■a 



- 0.2- 45 keV 
; EPICPN + M0S 



-41:10100 

0.2-45 keV 
EPIC PN + M0S 



IB 13 . 

' ' ■*» 

9 



.» 



~ 30 45O0 30 044)00 

ll-l l—l I 1—1 1—1 1—1 I 1—1 1—1 1—1 I 1—1 1—1 I I L 



47^00 ^ 30 ^ ^ 04^ OO ^ 30 



Figure 3. XMM-Newton EPIC M 3 1 images in the (0.2-4.5) keV XID band: North 3 (upper left), North 2 (upper right), North 1 (lower left), 
and Halo 4 (lower right). 



have originated if a too low background has been used for the 
long observations (exposures up to 50 ks) taken during the 
"reduced pointing phase" of the ROSAT mission. As Table 3 
shows this problem is not affecting the sources of the second 
ROSAT survey LH. There are many ROSAT PSPC sources de- 
tected with high LH that are missing in our XMM-Newton cat- 
alogue. Many of these sources have already been classified as 
transient comparing with the M 31 Einstein catalogue or just 
ROSAT survey I and II by SHP97 and SHL2001. At the posi- 
tion of one of these possible transients ([SHL2001] 240) we de- 
tect a source with a luminosity of ~1.5xl0 35 erg s _1 (if located 



in M 31), which is about a factor of 100 fainter than during the 
outburst detected by ROSAT. 

There are Chandra catalogues of M 31 X-ray sources 
based on observations of the centre area with the ACIS-I (204 
sources, 43 not detected by XMM-Newton, Kong et al. 2002b) 
and HRC (142 sources, 18 not detected by XMM-Newton, 
Kaaret 2002) detectors, and of shorter observations rastering 
part of the disk with the HRC (166 sources, 4 outside XMM- 
Newton field, 23 not detected by XMM-Newton, Williams et al. 
2004). While the Chandra centre area surveys are fully cov- 
ered by the XMM-Newton survey, the Chandra disk raster 
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Figure 3. (continued) XMM-Newton EPIC M 31 images in the (0.2-4.5) keV XID band: South 1 (left) and South 2 (right). The images are 
integrated in 2" pixels and smoothed with a Gaussian of FWHM 10", corrected for un-vignetted exposure and masked for exposures above 5 
ks for the individual cameras. Contour levels are (4, 8, 16, 32, 64, 128) x 10~ 6 ct s -1 pix~' including a factor of two smoothing. Sources from 
the catalogue are indicated. 



Figure 4. Logarithmically-scaled, three- 
colour XMM-Newton EPIC low back- 
ground image of the M 31 medium and 
thin filter observations combining PN and 
MOS1 and MOS2 cameras. Red, green and 
blue show respectively the (0.2-1.0) keV, 
(1.0-2.0) keV and (2.0-12.0) keV bands. 
The data in each energy band have been 
smoothed with a Gaussian of FWHM 20" 
to use an average point spread function 
for the different off-axis angles. The cam- 
era images have been corrected for un- 
vignetted exposure and masked for expo- 
sure above 5 ks. The image scale and the 
optical D25 ellipse of M 31 are marked. 
Shown in detail to the upper right is a fac- 
tor three zoom-in to the bulge area using 
images binned to 1" and smoothed with a 
Gaussian FWHM 5" corresponding to the 
point-spread-function (PSF) in the center of 
the field-of-view. 



partly exceeds the field covered by the XMM-Newton obser- 
vations. In addition to these larger Chandra survey papers, 
there are shorter lists of e.g. G1C sources and SSSs and re- 
ports of transients detected by Chandra. Correlations are given 
in the XID column of Table 2 and further references in the 
footnotes. While most of the Chandra detected sources can 
be resolved with XMM-Newton EPIC, there are sources close 
to the M 31 nucleus (319,321,325) or with arcsecond separa- 
tion (360,366) that remain unresolved. Many of the Chandra 



sources not detected by XMM-Newton are rather faint (lumi- 
nosities of a few 10 36 erg s _1 in the centre area and down to 
10 35 erg s _1 further out in the bulge) and therefore not resolved 
from the surrounding diffuse emission. The remaining unde- 
tected Chandra sources are rather bright and all but one of 
them (source nl-77 from Williams et al. (2004)) are already 
classified as transient (sl-79,sl-80,r3-46,sl-82,rl-34,r3-43,r2- 
28,rl-34,sl-85,rl-23,rl-19,nl-85). 
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Table 3. Summary of cross-correlation with the M 31 catalogue of the first (July 1991) and second (August 1992) ROSAT PSPC survey 
(SSHP97,SHL2001). ROSAT sources, not detected in our XMM-Newton analysis, are arranged according to detection likelihood (LH) levels. 



SSHP97 source numbers 



SHL2001 source numbers 



outside XMM-Newton field: 



resolved by XMM-Newton: 



not detected, LH<12: 



not detected, 12<LH<15: 



not detected, LH>15: 



1-24,27-33,35,37-41,46-48,50,52,54,58,60, 

61,65,70-72,78,79,81,85,87,88,91,95,98-100, 

102,104,110,114,120,122,123,125,129,134, 

135,137,148,154,159,169,171,180,183,185, 

1 99,202,209,2 1 2,2 1 6,2 1 9,224,227,230,236, 

237,242,245,246,257-260,264,265,268,269, 

271,273-276,282,285-287,290,293,294,298, 

300,305,312,314,318-320,327,335,339,341, 

342,350,356,358,362-365,367,371,373,374, 

376,377,381-383,385-396 

131,177,186,198,200,211,220,250,370 
(12 of 18 sources) 

59,62,63,69,84,86, 1 1 3, 149, 1 6 1 ,307,329,330 
(9 of 29 sources) 

49,82,93,109, 128, 196,283,334,372 
(22 of 193 sources) 

68(LH=14294),77(25),80(16),126(217), 

133(40),160(26),166(17),190(113),191(55), 

192(54),203(103),214(400),215(251), 

232(104),263(38),270(40),277(16),284(20), 

309(82),325(21),331(20),340(28) 



1-32,34,35,37^13,45,47,50,52,53,55,57,58,60-63, 

66,67,69,74,75,77,78,80,81,83,85,87-91,93,98,100, 

103,106,107,109,113,115-120,123,125,127,128,130, 

131,133-135,141,145,146,149,155,158,160,164,165, 

167,176,179,181,185,186,189,192,196,201,202,204, 

205,210,215,220,221,224,225,228,233,234,239,243, 

245,248,251-253,259,260,263-265,270,271,274,276, 

278,280-282,284,290-293,295-298,300,302,305,309, 

311,313,317,325,326,328,336,355,358,361,364-367, 

369,371-373,375,377-379,381-383,387-396 

129,177,178,183,191,193,195,197,206,214,254,272, 
374 

(3 of 21 sources) 
96,238,307 

(2 of 17 sources) 
76,79 

(18 of 167 sources) 

104(LH=901),121(94),171(43),173(317),190(216), 
207(98),208(299),230(76),232(1166),246(40), 
256(60),267(22),322(2703),324(148),344(41), 
356(15),380(17),384(16) 



Chandra spatially resolved several SNRs in M 31. Only 
some of them are detected by XMM-Newton (see Sect. 6.4). 

Only lists of bright XMM-Newton sources were published 
up to now (see e.g. Osborne et al. 2001; Trudolyubov et al. 
2004). All of these sources are contained in our catalogue. For 
bright sources in M 31, many photons were collected by the 
XMM-Newton EPIC cameras resulting in detailed light curves 
and spectra that could be analyzed in detail and allowed us to 
define the source class (see references in Table 2 and the dis- 
cussion of source classes in the following section). 

6. Classes of point-like X-ray sources detected in 
the direction of M 31 

To identify the X-ray sources in the M 31 field we searched 
for correlations around the X-ray source positions within a ra- 
dius of 3 x (cr s , a , + (Tsyst) in the SIMBAD and NED archives 
and within several catalogues. In columns 68 to 73 of Table 2, 
we give extraction information from the USNO-B1 catalogue 
(name, number of objects within search area, distance, B2, R2 
and I magnitude of the brightest object). To improve on the reli- 
ability of identifications we have used the B and R magnitudes 
to calculate log(<&-) = log(/ x ) + (m B2 + m R2 )/(2 x 2.5) + 5.37, 

r opi 

following Maccacaro et al. (1988, see Col. 74). 

The X-ray sources in the catalogue are identified or classi- 
fied based on properties in the X-ray (HRs, variability, extent) 
and of correlated objects in other wavelength regimes (Table 2, 
Cols. 75, 76). The criteria are summarized in Table 4 and are 
similar to the ones used in Paper I for M 33. However, one 



has to keep in mind, that M 31 is more massive than M 33 
with more interstellar matter in the disk which in addition is 
more inclined. Therefore sources in or behind M 3 1 may suf- 
fer higher absorption than similar sources in the M 33 field. 
To take this into account, we compare our HRs to model HRs 
determined assuming Galactic foreground absorption and ad- 
ditional absorption within M 31 of 9xl0 21 ctrT 2 , a factor of 2 
higher than in Paper I. We count a source as identified, if at 
least two criteria secure the identification. Otherwise, we only 
count a source as classified (indicated by pointed brackets). 

We plot X-ray colour/colour diagrams based on the HRs 
(see Fig. 5). Sources are plotted as dots if the error in both 
contributing HRs is below 0.2. Classified and identified sources 
are plotted as symbols even if the error in the contributing HRs 
is greater than 0.2. Symbols including a dot therefore mark the 
well defined HRs of a class. To identify areas of specific source 
classes in the plots, we over-plot colours of sources, derived 
from measured XMM-Newton spectra and model simulations. 

Identification and classification criteria and results are dis- 
cussed in detail in the subsections on individual source classes 
below. Many foreground stars, SSSs and SNRs can be classified 
or identified. However, besides a few clearly identified XRBs 
and AGN, and SNR candidates from positions in other wave- 
lengths, we have no clear hardness ratio criteria (see below) 
to select XRBs, Crab-like SNR or AGN. They are all "hard" 
sources and we therefore introduced a class <hard> for sources 
with HR2 minus HR2 error greater than -0.2. In the follow- 
ing subsections we first discuss foreground and background 
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Table 4. Summary of identifications and classifications. 



Source type Selection criteria identified classified 



fg Star lo g(^) < - L0 and HR1 - EHR2 < °- 3 and HR3 ~ EHR3 < ~ 0A or not defined 6 90 

AGN Radio source and not classification as SNR from HR2 or optical/radio 1 36 

Gal optical id with galaxy 1 

GC1 X-ray extent and/or spectrum 1 1 

SSS HRI < 0.0, HR2 - EHR2 < -0.99 or HR2 not defined, HR3, HR4 not defined 18 

SNR HRI > -0. 1 and HR2 < -0.2 and not a fg Star, or id with optical/radio SNR 21 23 

G1C optical id 27 10 

XRB optical id or X-ray variability 7 9 

hard HR2 - EHR2 > -0.2 or only HR3 and/or HR4 defined, and no other classification 567 



sources in the catalogue and then sources within M 3 1 . Thirty- 
eight sources remain unidentified or without classification. 

6.1. Foreground stars (fg Star) 

Foreground stars are a class of X-ray sources that is homoge- 
neously distributed over the FOV of M 3 1 . The good position- 
ing of XMM-Newton and the available catalogues USNO-B1 
and 2MASS allow us to effectively select this type of sources. 
All our foreground stars and candidates are 2MASS sources. 
We identified two sources (139 and 709) with stars of known 
type in the SIMBAD data base. Two sources (552 and 565) 
are W UMa type variables identified in the DIRECT program 
(Kaluzny et al. 1998; Stanek et al. 1999). Two sources (284 and 
385) were globular cluster candidates that have been identified 
as foreground stars (Galleti et al. 2004). Their X-ray fluxes 
are in the range expected for these sources. In addition, we 
classify 90 sources as foreground stars (log(-j&-) < -1, see 

Jopl 

Maccacaro et al. (1988) and in addition HRI - EHR2 < 0.3 
and HR3 - EHR3 < -0.4). For several of the star candidates 
we estimate the type from the optical colours in the USNO- 
Bl catalogue using the stellar spectral flux library from Pickles 
(1998). For many candidates we can not determine the type 
from the optical colours. This may indicate that these sources 
are not isolated stars but more complicated systems or even in 
some cases galaxies. This expectation can for instance be con- 
firmed for source 766, which is resolved into two sources in 
DSS2 images. 

Several of the foreground star candidates close to the cen- 
tre of M 31 (219,231,304,378,423,434) have no entry in the 
USNO-B1 catalogue, however, they are clearly visible on 
Digital Sky Survey images, they are 2MASS sources and fulfill 
the X-ray hardness ratio selection criteria. Therefore, we also 
classify them as foreground stars. 

For three sources (2,802,835) the HRs indicate a fore- 
ground star, however, log(-£-) are (-0.98,-0.87,-0.81) slightly 

/opt 

bigger than the assumed limit for classification. Comparing 
USNO-B1 and 2MASS optical magnitudes for these sources 
clearly indicates variability. Therefore, we still include them as 
foreground star candidates. 

Searching for correlations in the SIMBAD database, we 
found that source 136 correlates with the recurrent nova Rosino 



61. Rosino (1973) reports: "maxima of 16.4 pg and 17.7 pg 
have been respectively observed on August 12, 1966 and Oct 
25, 1968. In both cases, the star was near maximum for only 
one or two days, rapidly declining below the limit of visibility. 
Although the apparent magnitude at maximum may be consis- 
tent with that of a nova at the distance of M 3 1 , the light curve 
is abnormal, even for a recurrent nova. An alternative hypothe- 
sis is that the star may be not a nova, but a foreground U Gem 
variable, projected by chance over M 31." The X-ray hardness 
ratios of the source and its 2MASS classification support the 
latter explanation and we therefore add the source as a fore- 
ground star candidate. 

For seven sources (40,78,114,206,495,562,686), log(#-) 

/opt 

points at a stellar identification. However, the HRs including er- 
rors are outside the assumed limits for foreground stars. Source 
686 in addition correlates with a radio source making it a good 
candidate for an AGN behind M 3 1 . 

Optical spectroscopy of all foreground star candidates is 
needed to prove the suggested identification. 

6.2. Galaxies (Gal), galaxy clusters (GCI) and AGN 

Already after the Einstein observations the X-ray bright source 
[TF91] 51 (315) was identified with the local group galaxy 
M 32. Mainly based on ASCA and ROSAT HRI observations, 
Loewenstein et al. (1998) argued that the X-ray emission is 
dominated by a LMXB, slightly offset from the black hole nu- 
cleus of M 32. Chandra observations allowed to resolve this 
central emission of M 32 into three sources, one coincident 
with the nucleus and the by far brightest source at an offset 
from the galaxy nucleus of 8'.'3 coincident with the position 
of the proposed LMXB (Ho et al. 2003). In addition to these 
sources, Williams et al. (2004) detected a bright X-ray tran- 
sient in M 32 (sl-85). The XMM-Newton source 315 is posi- 
tioned far off-axis in observation si and does not allow us to 
resolve the three Chandra sources close to the nucleus. The 
luminosity of 5xl0 37 erg s" 1 and HRs of the source are cer- 
tainly dominated by the LMXB and typical for this source class 
(HRI = 0.31, HRI = 0.07, HR3 = -0.28, HR4 = -0.33). The 
M 32 transient detected by Chandra, was not active during the 
XMM-Newton observation. 
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Figure 5. Hardness ratios of sources detected by XMM-Newton EPIC. 
Shown as dots are only sources with HR errors smaller than 0.20 on 
both HR(i) and HR(i +1). Foreground stars and candidate are marked 
as big and small stars, AGN and candidates as big and small crosses, 
SSS candidates as triangles, SNR and candidates as big and small 
hexagons, GlCs and XRBs and candidates as big and small squares. 
In addition, we mark positions derived from measured XMM-Newton 
EPIC spectra and models for SSSs (S 1 to S4) as filled triangles, low 
mass XRBs (LI and L2) as filled squares, SNRs (N132D as Nl, 
IE 0102.2-7219 as N2, N157B as N3, Crab spectra as CI and C2) 
as filled hexagons, AGN (Al and A2) as asterisk. 



The brighter of the two extended sources in our survey 
(747) has been identified as a cluster of galaxies (GC1) at a 
redshift of z=0.293 based on the X-ray spectrum (Kotov et al. 
2003). The second extended source (832) is detected in the 
northeast part of the n3 field outside the optical D 2 s ellipse of 
M 31. We therefore classify this source as a GC1 candidate. 

There are no correlations with AGN with known redshift. 
However, the brightest X-ray source in the M 31 field (50) 
which was always active since the first Einstein observations 
correlates with a radio source and an unresolved optical object 
with a magnitude in the R band of 18.2. It has been identified 
as a BL Lac type AGN (NED). 

In addition, we classify 36 sources as AGN based on 
SIMBAD, NED, and other radio source correlations (NVSS, 
Braun (1990), Gelfand et al. (2004)) with the additional condi- 
tion of being not a SNR or SNR candidate from the X-ray HR. 
A final decision will only be possible from optical spectra as 
some of the sources may still be plerion type SNRs. In Fig. 5 
we include typical HRs for an AGN spectrum (power law with 
photon index of 1.7 assuming Galactic foreground absorption 
(Al) and absorption of 9xl0 21 ctrT 2 through M 31 (A2)). The 
AGN candidates populate the area in the HR diagrams expected 
from the model spectra. Many of the other sources in that HR 
diagram area - now just classified as hard - may turn out to be 
AGN. 

6.3. Super-Soft Sources (SSSs) 

SSSs show black body spectra with temperatures below 50 eV, 
radiate close to the Eddington luminosity of a 1 Mq object and 
are believed to be white dwarf systems steadily burning hy- 
drogen at the surface. They were identified as a class of X-ray 
sources by ROSAT and are often observed as transient X-ray 
sources (see Greiner 2000, and references therein). In the cat- 
alogue, SSSs are only classified using their HRs. To guide the 
classification we calculated SSS HRs assuming a 25 and 50 eV 
black body spectrum assuming Galactic foreground absorption 
(SI and S3) and absorption of 9xl0 21 ctrT 2 within M 31 (S2 
and S4). At the high absorption, the sources would no longer be 
detectable with XMM-Newton as the intrinsic flux (expected to 
be around 10 38 erg s~') would be reduced by 5 and 4 orders of 
magnitude for the two model spectra, respectively. This led to 
the selection criteria in Table 4 and the classification of 1 8 M 3 1 
SSSs. They are detected with absorbed XID band luminosities 
of 6xl0 34 erg s _1 to 9xl0 36 erg s" 1 . Only eight of the sources 
and the model values can be plotted in Fig. 5. For the other 10, 
HR2 is not determined. 

One of the SSS candidates (430) was already detected by 
Einstein and ROSAT. It was found to be variable, however 
not classified as a SSS in the ROSAT survey which may be 
caused by nearby sources unresolved by the ROSAT PSPC. 
This source as well as another six (313,320,336,359,369,431) 
were also detected by Chandra. Sources 320,336,367,431 were 
classified as supersoft transients (Williams et al. 2004; Di 
Stefano et al. 2004). Source 431 shows a 865s period in the 
XMM-Newton data reminiscent of a rotation period of a mag- 
netized white dwarf (Osborne et al. 2001). Four of the sources 
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correlate with optical novae (313,347,359,456) and two cor- 
relate with 2MASS sources, indicating a dusty surrounding 
(191,401). 

A special case is source 352, where HRl indicates a SSS, 
however, HRl is 0.94 and within the error not compatible with 
-1.0 and HR3 and HR4 are not undefined. Therefore we do 
not classify the source as a SSS. On the other hand, this source 
was classified as a variable SSS from Chandra observations 
(Di Stefano et al. 2004). It is detected by Einstein and ROSAT. 
There are several possibilities: Due to the position close to the 
centre of M 31 and the high source density in this region, 
XMM-Newton may not be able to resolve this source from 
nearby sources, the source detection may be fooled by the sur- 
rounding diffuse emission, or it is not a source with typical SSS 
properties. 

Another source (395), classified as a transient SSS from 
Chandra observations (Williams et al. 2004; Di Stefano et al. 
2004), has <hard> HRs in XMM-Newton and is classified here 
as XRB candidate. 

In the discussion above already two sources (352,395) were 
mentioned that are also covered by the Chandra survey for 
SSSs and QSSs. In this survey, Di Stefano et al. (2004) deter- 
mined 33 SSSs and QSSs using rates and hardness ratio crite- 
ria based on the Chandra ACIS-S (0. 1-1 . 1), ( 1 . 1-2), and (2-7) 
keV bands. As the authors discuss, this survey not only selects 
classical SSSs as defined above but also includes sources with 
effective temperatures as high as 350 eV. This expectation is 
confirmed by the 16 sources of the catalogue which coincide 
with sources from the XMM-Newton catalogue. Here, three 
of them are classified as SSS candidates (320,336,369), one 
as SNR (154), three as SNR candidates (295,318,704), four 
as foreground star candidates (128,157,231,727), two as XRB 
candidates (321,395), and one even as hard (169). 

6.4. Supernova remnants (SNR) 

SNRs can be separated in sources where thermal components 
dominate the X-ray spectrum below 2 keV (examples are 
N132D in the Large Magellanic Cloud and IE 0102.2-7219 
in the Small Magellanic Cloud) and so called "plerions" with 
power law spectra (examples are the Crab nebula and N157B 
in the Large Magellanic Cloud). To guide the classification we 
calculated HRs from archival XMM-Newton spectra of these 
SNRs. Spectra of N132D (Nl), IE 0102.2-7219 (N2), and 
N157B (N3) can directly be compared to M 31 SNRs as they 
are seen through comparable foreground absorption. For the 
Crab nebula spectrum, we assumed Galactic foreground ab- 
sorption (CI) and absorption of 9xl0 21 cirT 2 within of M 31 
(C2). It is clear from Fig. 5, that "thermal" SNRs are located 
in areas of the X-ray colour/colour diagrams that only overlap 
with foreground stars. If we assume that we have identified all 
foreground star candidates from the optical correlation and in- 
spection of the optical images, the remaining sources can be 
classified as SNRs with the criteria given in Table 4. Compared 
to our M 33 analysis in Paper I, we relaxed our HRl selection 
a bit to adjust to the higher absorption in the M 31 disk. 



Extensive searches for SNRs in M 31 were performed in 
the optical (Dodorico et al. 1980; Magnieret al. 1995, and ref- 
erences therein) and radio (Braun & Walterbos 1993). Before 
Chandra and XMM-Newton, M 3 1 X-ray sources were iden- 
tified as SNR just by positional coincidence with optical and 
radio catalogues. Also for our catalogue, many X-ray sources 
correlate with optical and/or radio SNRs or radio sources (see 
Table 2). With Chandra and XMM-Newton, we now have the 
possibility to identify M 3 1 SNRs by their X-ray properties (ex- 
tent, HRs, or spectra). Chandra observations of sources 260, 
342 and 454 (Kong et al. 2002a; Williams et al. 2004) were 
analyzed in detail and the SNR nature determined by X-ray 
extent and optical/radio correlation. We count sources as SNR 
identifications (21) if they correlate with optical and/or radio 
SNRs or radio sources and fulfill the HR criteria (see above). 
We count sources as classified SNRs (23) if they either fulfill 
the HR criteria (22) or correlate with an optical/radio SNR (1). 
The SNRs and candidates cover an absorbed XID band lumi- 
nosity range from 4xl0 34 erg s _1 to 5xl0 36 erg s _1 . The 22 HR 
selected SNR candidates have just been selected from their X- 
ray properties and significantly enlarge the number of M 31 
SNR candidates. For three of these sources, the classification 
may be questionable: From Chandra observations (Kong et al. 
2002a), source 316 is classified as variable, and 295 and 318 
correlate with a star of the catalogue by Haiman et al. (1994) . 

A special case is source 642 with an XID band luminos- 
ity of l.OxlO 36 erg s _1 , which correlates with the optical SNR 
candidate [MPV95] 1-013 (Magnieret al. 1995). The XMM- 
Newton HRs clearly point at a hard spectrum as expected for 
a plerion embedded in the M 31 disk. If the identification is 
correct this source would be the first plerion detected outside 
the Galaxy or the Magellanic Clouds. Final proof can only be 
achieved if the source turns out to be extended and shows no 
time variability. 

Kong et al. (2003b) reported the Chandra detection of two 
additional resolved X-ray SNRs in the centre of M 3 1 . The first 
of the sources lies about 40" northeast of the XMM-Newton 
source pair 338,341, the second inbetween these two sources. 
With X-ray luminosities of (4 and 8)xl0 35 erg s _1 in the (0.3- 
7) keV band, they are too faint to be detected by XMM-Newton 
in the neighbourhood of the bright sources and surrounding dif- 
fuse emission. 

6.5. Globular cluster sources (GIC) 

A significant part of the luminous X-ray sources in the Galaxy 
and M 31 are found in globular clusters. Most of the GlCs in 
the Galaxy show bursts and therefore are low mass neutron star 
X-ray binaries (see e.g. Verbunt & Lewin 2004). There have 
been extensive surveys for globular clusters in M 3 1 in the op- 
tical and infrared band (see Magnier 1993, Galleti et al. 2004). 
From the ROSAT survey sources, SHL2001 identify 33 sources 
as GlCs. We correlated the XMM-Newton catalogue with the 
catalogues of Galleti et al. (2004) and Magnier (1993) and add 
one candidate, found in SIMBAD. We count sources as identi- 
fications (27) when the correlating source is listed as confirmed 
GIC in Galleti et al. (2004). The remaining correlations are 
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counted as classified. Two sources (91,146), correlating with 
G1C candidates in Galleti et al. (2004) showed X-ray HRs typ- 
ical of SNRs, source 91 correlates with an optical SNR candi- 
date and both with radio sources. We therefore identify them as 
SNRs. 

In Fig. 5, the GlCs are plotted with the same symbols as 
the XRBs. All of our G1C identifications and candidates have 
been reported in earlier work and many of them have been clas- 
sified as time variable in X-rays (see Table 2). They cover an 
absorbed XID band luminosity range from 4.5xl0 35 erg s _1 to 
2.4x1 38 erg s _1 . Only the brightest source (351) has a lumi- 
nosity that seems to be at the upper end, allowed for neutron 
star LMXBs. 

For source 414, Trudolyubov et al. (2002b) report inten- 
sity dips with a period of 2.78 hr reminiscent of a neutron star 
LMXB. 

6.6. X-ray binaries (XRB) 

As already mentioned in the introduction to this section, ex- 
pected spectra of XRBs are similar to AGN and Crab-like 
SNRs. To guide the classification we calculated HRs as ex- 
pected for low mass XRBs (5 keV thermal Bremsstrahlung 
spectrum assuming Galactic foreground absorption (LI) and 
absorption of 9xl0 21 cirT 2 within M 31 (L2)). High mass XRB 
(HMXB) spectra are expected to be even harder. As can be seen 
in Fig. 5 the different source classes do not separate. 

Detailed work on individual sources using XMM-Newton 
and/or Chandra data, has identified four black hole XRBs 
(sources 257,287,310,384, Barnard et al. 2003b; Trudolyubov 
et al. 2001; Barnard et al. 2004), three neutron star LMXBs 
(297,353,403, Barnard et al. 2003b; Kaaret 2002; Mangano 
et al. 2004) and an XRB pulsar (544 Trudolyubov et al. 2004). 
In addition, we classify X-ray transients from XMM-Newton 
and/or Chandra data as XRB candidates. While in general 
transient behaviour of bright X-ray sources indicates an XRB 
nature, in this selection we may also pick up variable back- 
ground sources. Such a source may be (405) which also corre- 
lates with a radio source. The XRBs selected in this way cover 
the same area in the X-ray colour/colour diagrams as the GlCs 
(see Fig. 5). The absorbed XID band luminosities range from 
8.4xl0 35 erg s" 1 to 2.8xl0 38 erg s -1 . 

Besides the Chandra transients covered in the XMM- 
Newton catalogue, we identified many transients in earlier cat- 
alogues in Sect. 5 that are further XRB candidates. There cer- 
tainly are more transients XRBs buried in the XMM-Newton 
catalogue which only could be discovered by a detailed com- 
parison of fluxes of all the XMM-Newton sources with ear- 
lier missions. Such an analysis was outside the scope of this 
work and, due to the lower sensitivity of the catalogues be- 
fore XMM-Newton and Chandra, transient sources detected 
on this basis may be less reliable. On the other hand, also X- 
ray sources classified as variable may turn out to be XRBs. 
However, this criterion is not unique as it would also cover 
many variable background sources. Up to now no high mass 
XRBs are identified in M 3 1 . Several identifications with emis- 



sion line objects (EmO) in SIMBAD may be good candidates 
for Be-type high mass XRB. 

All the candidates mentioned above need further X-ray 
work to confirm an XRB nature. 

7. Summary 

We present a catalogue of 856 X-ray sources based on archival 
XMM-Newton observations covering an area of 1.24 square 
degree. We correlate our sources with earlier M 31 X-ray cat- 
alogues and use information from optical, infra-red and radio 
wavelengths and in addition their X-ray properties for source 
classification. We applied a source classification scheme simi- 
lar to the one successfully used in Paper I for M 33. As M 31 
sources we detect 21 SNRs and 23 SNR candidates, 18 SSS 
candidates, 7 XRBs and 9 XRB candidates, as well as 27 GlCs 
and 10 G1C candidates, which most likely are low mass XRBs 
within the G1C. Comparison to earlier X-ray surveys reveals 
transients not detected with XMM-Newton, which add to the 
number of M 31 XRBs. The number of 44 SNRs and candi- 
dates more than doubles the X-ray detected SNRs. 22 sources 
are new SNR candidates in M 3 1 based on X-ray selection cri- 
teria. Another SNR candidate may be the first plerion detected 
outside the Galaxy and the Magellanic Clouds. On the other 
hand, six sources are foreground stars and 90 foreground star 
candidates, one is a BL Lac type AGN and 36 are AGN can- 
didates, one source coincides with the Local Group galaxy M 
32, one with a background galaxy cluster and another is a GC1 
candidate, all sources not connected to M 3 1 . 

There are 567 sources classified as hard, which may either 
be XRBs or Crab-like SNRs in M 31 or background AGN. 
If source variability is detected SNRs may be excluded from 
source identification. Thirty-eight sources remain unidentified 
or without classification. 

The archival XMM-Newton M 3 1 observations allowed us 
to probe the point source population of the covered area signifi- 
cantly deeper than during the ROS AT PSPC surveys. However, 
the XMM-Newton EPIC coverage is rather inhomogeneous 
and a significant part of the M 3 1 disks has not been observed 
at all. To get a full census of the M 3 1 X-ray point sources it 
would be very important to fully cover this nearby galaxy with 
XMM-Newton observations. 
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